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Abstract 

An addition rule of impure density operators, which provides a pure state density operator, is for- 
mulated. Quantum interference including visibility property is discussed in the context of the density 
operator formalism. A measure of entanglement is then introduced as the norm of the matrix equal to 
the difference between a bipartite density matrix and the tensor product of partial traces. Entanglement 
for arbitrary quantum observables for multipartite systems is discussed. Star-product kernels are used 
to map the formulation of the addition rule of density operators onto the addition rule of symbols of the 
operators. Entanglement and nonlocalization of the pure state projector and allied operators are dis- 
cussed. Tomographic and Weyl symbols (tomograms and Wigner functions) are considered as examples. 
The squeezed-states and some spin-states (two qubits) are studied to illustrate the formalism. 

keywords: purification, entanglement, star-product, superposition 

1 Introduction 

Superposition principle of quantum states plays a key role in such physical phenomena as interference of 
matter waves Wave properties of electron are connected with de Broglie wave length expressed in terms 
of particle momentum . These properties are naturally described by a wave function associated with the 
particle's quantum state and obeying Schrodinger equation ||^. For a system with several degrees of freedom, 
the possibility to consider two subsystems — first one connected with some of degrees of freedom and the 
second one with the rest degrees of freedom, respectively, the superposition principle provides a construction 
of entangled states Discussing two subsystems of a given system implies the physical possibility to 
measure characteristic properties distinguishing the subsystems. 

Entangled states are the states which are constructed as a superposition of states, each of which has 
the wavefunction expressed as a product of wavefunctions depending on the different degrees of freedom. 
The mixed states of quantum systems are described by density operator Q. The superposition principle 
of pure quantum states has been formulated in [6-^8] in terms of a new addition rule of the density opera- 
tors. This addition rule corresponds to a purification procedure of a mixed quantum state obtained by the 
standard addition rule of the density operators. Relation of the purification procedure to reconstructing the 
entanglement structure of the mixed state of a bipartite system has been preliminarily discussed in . 

Various notions of measure of entanglement were suggested in [9-13]. All these measures are related to 
some operators associated to a bipartite quantum system. The aim of our work is to give the new addition 



rule of density operators describing the superposition of impure density matrices and this analysis generalizes 
the results of [6-8] where the coherent addition rule of pure density operators was formulated. We also define 
the measure of entanglement of bipartite and multipartite quantum systems considering intrinsic properties 
of the density operator describing the state. 

Since the density operators can be considered using different representations for their symbols like, for 
example, Wigner function [Q , Husimi-Kano function jl^, |l^ as well as singular quasidistribution jl^, |l^ , 
we discuss the addition rule of the density operators in terms of the addition rule of their symbols. To 
do this, we discuss the star-product of the operator symbols (see, for example, [19-23]). We consider also 
the case of density operator representation by the standard tomographic probability distribution, which is 
used to give a "probability" formulation of quantum mechanics [Q . We also give the formulation of both 
aspects of superposition principle, namely, the coherent addition of impure density operators and measure 
of entanglement in terms of the star-product quantization procedure. 

The paper is organized as follows. 

In section 2 the basic ideas of the construction of purification procedure and measure of entanglement 
are described. In section 3 the review of addition rule of pure density operators is presented including 
the presence of visibility parameter. In section 4 a new formula for purification of sum of impure density 
operators is obtained. In section 5 a short review of the star-product formalism is given. The purification 
formula for a symbol for an arbitrary kind of density operator which is obtained by purification of sum of the 
symbols of impure density operators is obtained in section 6. The example of the Wigner function addition 
rule in terms of star-product kernel is presented in section 7. Addition rule of tomographic symbols is 
considered in section 8. In section 9 the notions of entanglement and measure of entanglement are discussed 
in terms of intrinsic properties of the density operator of a composite system, while in section 10 a notion 
of entanglement is introduced for other observables. The example of two qubits is considered in section 11. 
A measure of entanglement of multimode squeezed state is presented in section 12. Purification procedure 
for separable density matrix is discussed in section 13. In section 14 the role of fiducial projector used to 
formulate the superposition principle in terms of density operators is considered. In conclusions (section 15) 
some perspectives are discussed, while in Appendix we give the proof of the theorem that for pure state pab 
of bipartite system AB eigenvalues and ranks of reduced density operators pA and pB are equal. 

2 General ideas 

The notion of dynamical variables generating an algebra of observables and the notion of states, which 
are dual to this algebra, are common to both classical dynamics and quantum dynamics; but in quantum 
dynamics the operators generating transformations form a noncommutative vector space while the classical 
algebra of dynamical variables is commutative. Since every true representation (realization) of a commutative 
algebra is one-dimensional, this is no longer true of a noncommutative algebra. This has the immediate 
consequence that while one can have states in which all dynamical variables have unique values in classical 
dynamics, this is not so in quantum dynamics. We have fundamental states in quantum dynamics in 
which all but a (commuting) set of variables have definite (dispersion-free) values for which other dynamical 
variables have a distribution of values. In the case where the eigenstates are continuous, this finds expression 
in generalized uncertainty relations [25-28]. The generalised uncertainty relations exist also for discrete 
observables like spin (see, e.g., p^). 

Classical states may be associated with distributions in phase space (that is, they have distributions 
in values for all dynamical variables) of which the extremal (pure) states are points in phase space (their 
distributions correspond to Dirac delta functions or Kronecker delta matrices for the fundamental dynamical 
variables). In contrast, the quantum states are linear functionals on the dynamical variables which map 
nonnegative operators into nonnegative numbers. A linear distribution that can be expressed as a convex 
combination (sum or integral) of other appropriate ones can be made up of extremal linear functionals which 
cannot be so decomposed. These extremal states are the 'pure states' of a quantum system. 

If all quantum states had distributions for some set of dynamical variables, how can we combine them 
to get dispersion free states for these dynamical variables? This depends on the characteristic quantum 
property of superposition of states. Thus, for example, in interference we compose two pure states to form a 
new pure state. In terms of states considered as linear functionals, the natural process is to generate convex 
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linear sums which are not pure states. But processes hke interference, diffraction and composition of Ught 
polarizations need a new procedure. 

An alternate formalism suitable for incorporating this composition law is given by the vector (Hilbert) 
space formalism — with each pure state we associate a vector (and its dual vector), and the expectation 
values are obtained by the dynamical variables acting on the state as linear operators, and forming the scalar 
product with the original dual. Real classical variables have their counterpart in selfadjoint linear operators 
and their expectation value in any state is real. Moreover, positive operators (nonnegative operators) have 
nonnegative expectation values. Since vectors permit linear combinations (over the field of complex numbers), 
superposition of states and interference are naturally explained. Since the sesquilinear tensor product of the 
vector and its dual is a linear functional, this can be compared with the formalism in terms of linear 
functionals. But the only linear functionals obtained as a sesquilinear form (outer product of a vector and 
its dual) are extremal. To get full correspondence with the linear functional formulation, we should form 
convex combinations of such outer products. Density matrices (nonextremal states) are convex combinations 
of bilinears in the vectors. 

Since the sesquilinear form (outer product) of states constitute pure states, we see that every 'state vector' 
corresponds to an extremal linear functional, which may be treated as a linear operator in the vector space, 
we should identify mixed states as nonextremal states of the convex set of states. But the correspondence of a 
state vector and its pure density operator is not one-to-one. The von Neumann 'measurement' is a projector 
to an eigenstate of particular self-adjoint dynamical variable. (In the case spectrum is continuous, von 
Neumann prescribes a nested family of spectral projections.) The prescription may be viewed as imparting 
a measure on the spectrum of the operator representing the dynamical variable. 

A generalization of this protocol is called POVM (positive operator valued measure). A von Neumann 
measurement results in a pure state of the system (or more generally a density matrix in the eigenspace 
of the operator). But a POVM may result in a mixed density matrix. The entire set of vectors {ipe^^} 
corresponds to a unique extremal state 

p = IpTp'^ . 

This extremal operator is idempotent, selfadjoint (of trace class) and satisfies 

= ip^'^^ip'^ = _ p 

since ^^tp = 1 for normalization. This equivalence class {e**'0} is identified as a ray — so extremal density 
matrices correspond one-to-one to rays in the vector (Hilbert) space. Rays do not constitute a vector space. 
Considered as a linear operator, since = p and Trp = T/jtV = 1, p is a projector of rank one. So extremal 
states are projectors associated with unique rays. A mixed state is a probabilitistic (nonnegative normalized 
linear) sum of projectors corresponding to sum of definite rays with definite probabilities, which may be 
chosen to be mutually orthogonal. 

In forming superposition of two vectors, their relative phases are important. So the rays by themselves 
are insufficient. They suffice to form mixtures. The question arises as to how to work with the projectors 
and yet get a superposition with definite (relative) phases within the formalism of density operators, either 
pure or mixed. 

This is accomplished using the method of 'purification' of a nonextremal state. Since there are many 
possible superpositions, we must have purification of an impure density to any of those states. This is 
accomplished by the use of a suitable fiducial projector; the choice of this projector determines the phases 
(or rather phase differences). We show that there are choices of the fiducial projector that can give any of 
the superpositions and that, even with one such superposition, the fiducial projector can be chosen amongst 
a continuous set of projectors. 

Having accomplished 'purification' we turn to another characteristic property of quantum system, namely, 
'quantum entanglement'. If we have a composite system, say AB^ composed of two subsystems A and 5, the 
generic states of AB contain information that is not obtained by considering the states of the subsystems. 
These may be refered to as nonlocal correlations between the subsystems A and _B, which cannot be attributed 
to causal connections. (We have processes which have significant relationship that cannot be accounted 
causally; somewhat like the notion 'synchronicity' by Carl Jung [^.) Schrodinger pointed out this as a 
characteristic property representing a quantum system. 

We can have correlations between dynamical variables measured in subsystems in classical dynamics. In 
a pure classical state, this automatically gives pure states with definite values for the subsystem variables. 
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But in the quantum system the situation is entirely different. A general pure state of the composite system 
gives impure states of the subsystems. In this case, the systems A and B are 'entangled'. For example, the 
singlet state of two spin-1/2 (particles) is entangled. Any spin component of the combined system AB gives 
a zero expectation value (singlet is total spin-0), the individual particles are completely unpolarized — any 
spin component has zero expectation value. Yet there is a definite correlation — if one spin is 'down', the 
other one is 'up' (and vice versa) with respect to any direction. 

We recognize that if an entangled pure state is considered as a state of two subsystems A and B, they 
could be mixed with the same rank but with corresponding eigenvalues and eigenprojectors. So the problem 
of recovering the original pure state from the two mixed (impure) states of the subsystems involves the 
restoration of nonlocal phase relations, characteristics of entangled states. We have evolved a method of 
restoring the pure state using an entangled fiducial projector. However, in this case, the restoration can be 
done if the subsystems have density matrices of the same rank and same eigenvalues. 

We may also recognize that operations on an entangled pure state may lead to an unentangled (Kronecker 
product) state for the composite system which remains pure when restricted to either subsystem. For 
example, if we act on the singlet state with the operator of the difference of two spins, we can obtain a 
triplet state which may be unentangled. The point is that the operator acting as the difference of the spins 
is itself 'nonlocal' in that it acts on both subsystems together. Similarly, we can take an unentangled pure 
state like the ±1 states for the total spin-1, and an entangled pure state can be obtained by either acting 
with some component of the difference of spins on the singlet state, or by a spin operator antisymmetric in 
the two spins on the entangled triplet state. Here again the operator has the ability to entangle. We shall 
make use of this property in generating entanglement. The entanglement and the relative phases are implicit 
in the fiducial projector chosen. 



3 Pure state addition 

In [6-8] the rule to add two pure states determined by their density operators pi and p2, which are projectors, 
was formulated. This rule corresponds to superposition of the state vectors ■01 and ■02 ■ This superposition 
of vectors is a standard tool to describe the quantum phenomenon of interference. The interference pattern 
is sensitive to relative phase of the two vectors. To describe the relative phase, a fiducial projector Pq was 
used. By means of the fiducial projector, the quantum interference can be described in terms of operators 
only without using the state vectors. The addition law of two orthogonal states reads (see [6-8] where the 
addition formula was written in slightly different form) 

P^PlPl+P2P2 + ====- (1) 

^yTT [P1P0P2PQ) 

where positive probabilities (numbers pi and P2) satisfy the normalization relation 

Pi+P2 = l. (2) 

One can check that the density operator (|^) is a projector, i.e. 

P^ = P P^ = P Trp = 1 p^ppi = pipi {i = 1, 2) (3) 

if Pi, P2 and Pq are projectors. 

We consider the case pip2 = and PqPi ^ and P0P2 7^ 0. 

The composition law (|l|) can be interpreted as the purification of the impure density operator 

Pim = Pipi + P2P2- (4) 

Nevertheless, the real meaning of the "purification" consists of the statement that for two given orthogonal 
projectors and given fiducial one which is not orthogonal to the given projectors, sum (|l|) is again a projector. 
Relation (|l|) can be extended to describe the purification procedure for the impure density operator of the 
form 

^^PkPk (5) 



k 
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where projectors are such that pkPm = (fc 7^ m) and positive probabihties pk satisfy the normahzation 
condition 

k 

Summation in (j^) and can be considered as summation over finite set of indices k or over infinite one for 
the system with infinite number of states. The generalization of (|l|) provides the purified density operator 

Formula (|l|) can be also extended to take into account that there is a visibility parameter 7 , a characteristic 
of the interference pattern with < 7 < 1. Equation (0) is generalized to the form 

P=PlPl+P2P2+l j====^. (8) 

For 7=1, equation (H) reduces to equation (|l|). 
For 7 = 0, one has the impure state (^). 

Thus, the visibility parameter 7 is a characteristic of completeness of the purification procedure of the 
density operator or of degree of decoherence of the initial pure superposition state (|l|). In the case of (|]) or 
(^, we have puppk = PkPk (no sum on k). 

Let us consider now two density operators of quantum states which are pure nonorthogonal states. In 
this case, one has normalization constant and the purification formula is 



P - 



, y/PlP2 (P1P0P2 +h-C.) 
PlPl + P2P2 -' 



i/Tr {P1P0P2P0 
The normalization constant J\f reads 



J\f-\ (9) 



^^^^^ ^ 2^p^e{TT{p,Pop2)) 
v/Tr ipiPoP2Po) 

4 Impure state addition 

We discuss now the addition rule of two mixed states. 

In the case of two impure states pi and p2, their sum can be decomposed in terms of orthogonal projectors 
Rn (i.e. satisfying i?„i?m = SnmRn)- 

P1PI+P2P2 ^^(^nRn ^W„ = l > 0. (11) 

n n 

One can consider formula ( pT|) as giving the result of mixture of pure states i?„ in (|^). Thus the procedure 
of addition of impure states can be fulfilled as follows. First, one writes the sum of impure states as a convex 
sum of orthogonal projectors and then carry out the purification given in equation (|^). 

In this case, one gets the pure density operator for the result of "deformed" addition of two impure 
states pi and p2 with probabilities pi and p2 given by the following expression which is a generalization of 
equation (pT]): 

RkPaRj /ir>\ 

where Rk are orthogonal eigenprojectors of the density operator and tUk are nonnegative eigenvalues of the 
density operator, i.e. 

(piPi + P2P2) Rk = (^kRk- (13) 
There is no sum on k. The fiducial projector Pq is chosen to satisfy the condition PoRk 7^ 0. 
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If one has addition of more than two impure states, i.e. the density operator of impure state has the 
form X^fLiPsPs, the result of 'deformed' addition rule has the same form as equation (p^, namely, 

where the eigenprojectors Rk and nonnegative eigenvalues ujk satisfy the equation 

(15) 



/ N \ 

( ^^PsPs J Rk = i^kRk (no sum on k). 



Thus the purification procedure which is expressed by the deformed addition rule denoted by the sign 
® in the left-hand side of equations (|l2|), (^4|) is reduced to obtaining eigenvectors (eigenprojectors) and 
eigenvalues of the nonnegative density operator ^^PsPs and applying ansatz with the fiducial projector 
Pq {PoRk 0) to construct the nonlinear expression in the right-hand side of equation (p^, (|^. This 
expression provides the purified density operator. 

In the limit case where the initial density operators ps are orthogonal projectors, equation ( p^ is reduced 
to equation (^ with obvious replacement ujk Pk- One can point out that only one projector Pq is sufficient 
to provide {N — 1) independent phase parameters in the case of addition of N orthogonal pure states. This 
projector must have nonzero overlap with all added pure state projectors. In the case of addition of N impure 
states, the number of independent phase parameters, which are contained in only one fiducial projector Pq, 
depends on the rank of the density operator PsPs and is one less than the rank. For clarity, we point 
out that the obtained addition formula for impure density operators assumes the existence of the fiducial 
projector. This existence is obvious from a geometrical point of view, but the explicit finding of this fiducial 
projector for given density operators is a different problem. 



5 Symbols and their star-product 

In quantum mechanics, observables are described by linear operators acting on the Hilbert space of states. 
In order to consider observables as functions, we review first a general construction and provide general 
relations and properties of a map from operators onto functions having in mind a map of density operator onto 
a function like Wigner distribution. Given a Hilbert space H and an operator A acting on this space, let us 
suppose that we have a set of operators J7(x) acting on H , where a rt-dimensional vector x = (xi, X2, ■ ■ ■ ,Xn) 
labels the particular operator in the set. We construct the c-number function /^(x) (we call it the symbol 
of operator A ) using the definition 

/^(x) = Tr [AuiK)] . (16) 
Let us suppose that there must exist a set of operators -D(x) such that 



A = / /^(x)i?(x) dx. 



(17) 



We will consider relations ( p^ and ( pT[ ) as relations determining the invertible map from the operator A 
onto function /^(x). 

The product (star-product) of two functions /^(x) and /^^(x) corresponding to two operators A and B 
is defined by the relations 



/Ib(^) - /iW * fsi^) Tr ABU{^) 



(18) 



Since the standard product of operators on a Hilbert space is an associative product, formula ( |18[ ) defines 
an associative product for the functions /^(x). 
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6 Superposition rule in terms of symbols of density operators 



Using formulas (|l6| ) and (|l7|), one can write down a composition rule for two symbols /^(x) and /^(x), 
which determines the star-product of these symbols bilinear in the two symbols, 



/^(x)*/^(x) = y'/4(x")/^(x')i^(x",x',x)dx'dx". 



(19) 



The kernel in the integral of ( p^ ) is determined by the trace of product of the basic operators, which we use 
to construct the map 



i^(x",x',x) ==Tr D{ji.")D(^')U(x) 



(20) 



Formula (^0|) can be extended for the case of the star- product of N symbols of operators Ai, A2, . . . , An . 
Thus one has 



(21) 



/ai W * /a^W * • • • * /a„ W = J /ii(xi)/A,(x2) • • • Ja^M 
xK (xi, X2, . . . , xjv, x) dxi dx.2 ■ ■ ■ dx^v 

where the kernel has the form 

if (xi,X2,...,xjv,x) =Tr[^(xi)i)(x2)---7)(xjv)?7(x) . 
The trace of an operator A^ is determined by the kernel as follows 

Tri^ = / /^(xi)/^(x2)---/^(x^) 



(22) 



xTr 



I?(X1)^(X2)---^(XA,) 



dxi dx2 • • • dx 



N 



Tr (AB) 



/^(xi)/^(x2)Tr ^(xi)i3(x2) 



dxi (iX2 . 



(23) 
(24) 



When the operator A is a density operator of a quantum state, formula ( p3| ) for N = 2 determines the purity 
parameter of the state. 

Formulas ( ^l[ ) and ( p2| ) can be used to formulate the addition law of the density operators of orthogonal 
pure states pi as addition law for their symbols fp.{x.). In the case of purification of the sum ^j^PkPk of 
pure states pk by means of a fiducial projector Pq, one has the following symbols: 

/p(x) for purified density operator, 

fp^ (x) for pure state with density operator pk , 

/o(x) for fiducial projector Pq. 

The formula describing quantum interference in terms of symbols of the density operators reads 



J,,-, , //pJxi)/o(x2)/p^.(x3)ii:(xi,X2,X3,x)dxidx2dx3 

' I fpk ) /o (X2 ) fp, (X3 ) /o (X4 ) fc (xi , X2 , X3 , X4 ) dxi dx2 dx3 dX4 



jk 



(25) 



where the kernel iir(xi, X2, X3, x) is defined by equation ( |22| ) while the kernel A;(xi,X2,X3,X4), which deter- 
mines trace of the product of N operators {N = 4) in terms of their symbols, reads 



fc(xi,X2, . . . ,XAr) = Tr l)(xi) L'(X2) . ..D(xn) 



(26) 



In the case of purification of the sum of impure states ([141), one has analogous formula with the replacement 
Pk Rk and pk^ ujk- 
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7 Weyl symbol addition and interference in terms of Wigner— 
Moyal functions 

In this section, we will consider a known example of Wigner distribution which is related to the Heisenberg- 
Weyl-group representation. The interference can be described in terms of Wigner functions using for- 
mula (25). The displacement operator 



D{ay^) — exp (cty^a) — ot^a) 

where 

ay^ — xi+ 1X2 a* = a;i — 1x2 
and real numbers xi and X2 are expressed in terms of position and momentum as 



xi = 



x/2 



X2 = 



P 

V2 



(27) 
(28) 

(29) 



determines the basic operators defining the Weyl map. Thus, one has for the basic operators of the map the 
following expressions [E2[: 



(30) 
(31) 



The operator (—1)" ° is the parity operator (—1)° " = P, with the matrix elements given in the position (or 
momentum) representation by the formula 

{x\P\y)^S{x + y). (32) 

The Weyl symbol of density operator p is defined by (|6|) where we use the operator (^^ and make the 
replacement / W, which is the state Wigner function, and it reads 

Wp(Q!) = 2Tr pD{a){-lf^D{-a) {a = a^). 

To describe the star-product of Weyl symbols, we introduce a generalization of notation (|29| ) 



— \-^ki 1 •^k2 ) ^ki — "^7^ 



Xk2 



Pk 

V2- 



Then 



ak = {qk + ipk) fc = 1, 2, . . . , TV. 
The kernel of the star-product of {N — 1) Weyl symbols has the form 

K (ai, a2, . . . , ajv) = Tr [j7(xjv)nfj/^(xfc) 



The kernel can be rewritten in terms of the complex numbers ai (i = 1, 2, . . . , N) as |22 
(ai, a2, . . . , QfTv) = 



N-1 JV-1 

j>i 1=1 



Af-l 



-N, 



(33) 



(34) 



where q = —1. 
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Kernel for trace of the product of four operators reads 



fc(ai, Q!2, as, 04) = { 



ai — a2 + as — a4) 



X exp 



I - 2 (aia 



2 — aia3 + 010:4 + 020:3 — 020^4 + 0304 j — c.c. 



(35) 



Having the above kernels we can obtain Weyl symbol of pure density operator (which we found by means of 
purification of mixture of several states) by inserting the kernels and the Wigner functions into (|2^). The 
explicit result for addition of two Wigner functions was given in B in a different form. 



8 Symplectic tomograms and superposition principle 

Now we consider the example of tomograms. 

According to the general scheme, one can introduce for an operator A the tomographic symbol /^(x), 
where x = [xi, X2, x^) = (X, /i,i/), which we denote here as w^(X, /i, j/) depending on the position X and 
the parameters fi and v of the reference frame |24| 



wxiX,n,iy) = Tr ij7(x) 



The operator U (x) is given by 



C/(x) = U{X, n,y) = 5{X ~ iiq- vp) = \X\5 1 



X X 



where q and p are position and momentum operators. 
The inverse transform will be of the form 



A 



w^{X, n, u)D{X, /i, v) dX djjL dv 



where 



Zj(x) = b{X, fi.v) ~ — exp {iX — ivp — ifjLq) . 
27r 



The kernel defining the star-product of two tomograms 

K{x!', x', x) = Tr \d{X", fi" , v")b{X' , fi' , v')U{X, ^, v) 

reads 



K{Xi,^l, l'l,X2, Ai2, J^2, X^, v) 

X exp 



47r2 



({Vifl2 - V2f^l) + + 2X2 - ^ ^ X 

2 \ V 

The trace of product of four basic operators, which provides kernel to caclulate the denominator in addition 
formula of density operators (ES) reads 



fc(xi,X2, . . . ,XAr) = Tr 



N 



nf.i^ (X,., iik. vu) = (27r)i-^<5 E M E 



Vk 



X exp 



I \fc=l fe<s = l 



\k=l 



\k=l 



VklJ-s - tJ-kVs) 



iV = 4. 



(36) 



Having the above kernels we can obtain the tomogram of pure density operator inserting the kernels into ( p5| ) 
and making there the replacement f w. In a different form the explicit result for addition of tomograms 
was done in |6(|. 
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9 Notion of entanglement 



Another quantum-mechanical property related to superposition principle of states in bipartite and multi- 
partite systems is entanglement. Let us have density operator pab of composite system AB which has two 
subsystems A and B. This means that there exist experimental possibilities to measure the properties of 
the subsystem A and of the subsystem B. The density operator pab determines two density operators of 
the subsystems 

PA = TrspAB and ps = Ttapab- 
Let us consider tensor product of the two subsystem density operators 

Paxb ^ Pa® Pb- 

There is difference of two density operators 

Rab = PAB - PA® Pb- 

This difference is a characteristic of entanglement. If the system is in the state pab, which is disentangled, 
the operator Rab = 0. Numerical characteristic of entanglement is described by nonzero matrix elements of 
the operator Rab- A basic independent (invariant) characteristic of the operator Rab is the number 

e = Tr{RABR\B)- (37) 
This number can be considered as a measure of entanglement. Since R\g — Rab one has 

e = Tr {r:\s) - 

There are other numerical characteristics of entanglement like traces of higher powers of the matrix Rab 

e(")=Tr(i?:^+i). 
The state pab is characterized by purity parameter 

PAB = Trp^s 
and the state paxb has its own purity parameter 

PAxB — PAPB 

where 

PA = T^t^Pa and pB = Trp|. 

Since 

Rab = Pab + (pa ® PbY - PabPa ® Pb - Pa® PbPab (38) 

one has 

Tri?iB - Trpis + Tr (p^ ® pi) - 2Tr [pab PA ® Pb)- (39) 
Thus we get the measure of entanglement in the form 

e = Pab + PaPb - 2^pabPaPb cos 9. (40) 

The last term in the right-hand side of ( |40| ) is determined as 

y/PABPAPB COS 6* = Tr {pabPa ^ Pb) - 

For pure state pab = 1, one has pA = Pb — P (see Appendix) and it gives e = \ + p^ — 2pcos6. The 
introduced angle 6 and parameters pA, Pb and pab can be functionally dependent. In view of this, finding 
maxima or minima of the entanglement measure needs taking into account this dependence. 
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In fact, the measure of entanglement ( p7|) is defined using the notion of distance between two density 
operators (see, e.g. |]30| ). The connection of distance with measure of entanglement is natural and it 
was discussed , e.g., in [31-34]. In the present work, we use the Hilbert-Scmidt distance as measure of 
entanglement but the novelty of suggested measure of entanglement is related to the choice of the density 
operators being compared. We use the distance between the system density operator and the tensor product 
of the partial traces over the subsystem degrees of freedom. This characteristic is intrinsic because it is 
contained in the state density operator only. The geometrical sense of the notion of the entanglement 
measure can be clarified using an analogy with distance between the points on Eucledean vectors |a — b|, 
where real vectors a and b describe the points. With this definition, a partially separabble system has a 
nonzero entanglement. 

Each matrix can be considered as a complex vector. The standard scalar product of any two vectors C 
and D can be always treated as 

s 

If one considers the matrix element of the two matrices C and D as components of the vectors C and D, 
one also has 

C • D = TrpC^). 
For Hermitian matrices C — and D — D\ 

C - D = Tr(Ci:i). 

If one considers a density matrix as the vector, purity parameter plays the role of square of the vector length, 
so one has for the purity parameter the inequality 

o < /i < 1. 

Description of the matrix Rab as a vector makes obvious that measure of entanglement (|3^ ) coincides with 
square of the vector length, which in turn is difference of two other vectors. This means that the length of 
the vectors under consideration which correspond to normalized density matrices is less than unity. Thus, 
the geometrical interpretation of measure of entanglement means that angle 9 in (|4^) is the angle between 
the two vectors. This angle can depend on the length of the vectors determining the purity parameters of 
the system and subsystems. The angle parameter is introduced only in order to illustrate the geometrical 
picture of the entanglement measure under discussion. 

10 Entanglement for arbitrary observables 

Usually the notion of entanglement is applied for the density operator. Mathematically the construction of 
measure of entanglement e given by ( |37| ) can be extended for arbitrary observable O represented in the form 
of sum of projectors 

O = ^ a„P„ 

n 

where the Pn are eigenprojectors and a„ are eigenvalues of the observable O, i.e. 

For a density operator, the eigenvalues are nonnegative numbers. For an arbitrary observable, the eigenvalues 
are real numbers and they can take negative values. If one has prescribed division of the system in terms of 
two subsystems A and B, the observable O which acts in Hilbert space of the system can be treated in the 
same manner as density operator in previous section. Thus one can define the reduced observables 

6 A = Trs 6 6b ^ Tr^ 6. 

The tensor product of the observables 

Oaxs ^6 a® 6b 
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acts in the Hilbert space of the system. 

The correlations of two subsystems captured by the observable O can be connected with a measure of 
entanglement like in the case of density operator. We define the measure of entanglement for the observable 
O as the number 



eo = Tr 



{O-Oa® Ob? 



Tr 



This number gives invariant description of a 'distance' between two Hermitian operators O and Oa ® Ob 
exactly in the same manner as in the case of distance between two density operators. 
Analogously, one can introduce positive parameters 



2 

kA 



where a^A and a^B are eigenvalues of the Hermitian matrices Oa and Ob, respectively. So formula ( ^ ) can 
be extended for arbitrary observable in the form 

eo = Mo + fJ'OAf^oB - 2^7*0/^0^7*08 cos 6* 
where we define cos 9 using the same geometrical interpretation of scalar product of two vectors 

^^oMoAMoB cos 6* = Tr ( 66 a <^ 6b 



Thus we introduced notion of entanglement for other Hermitian observables than for density operators. 
Of course, the inequalities for purity parameters in the case of density operators are not valid for other 
observables. 

One can make generalization introducing the measure of entanglement of fcth order of arbitrary observable 
to multipartite system AB . . . C using the definition of measure 



= Tr 



[O-Oj 



'Of 



'Ocf 



Oa = TrB...cO, 



/c = 2,3, 



.N. 



(41) 



For even k, the above parameter is nonnegative number. The measure can be normalized using the factor 



11 Example of two qubits 



Let us consider a density matrix with unit trace for two spins in the basis |t) and |J,) for the first spin and 
for the second spin, correspondingly, i.e., in the basis in four-dimensional space 

ITT) ITT) ITT) ITT). 



The Hermitian density matrix has a form 



The density matrix pA = Tr^p reads 

PA -~ 

and the density matrix pB — Tr^p reads 

PB = 



Pll 


Pl2 


Pl3 


Pl4 \ 


P21 


P22 


P23 


P24 


P31 


P32 


P33 


P34 


Pil 


PA2 


P43 


P44 / 


Pll 


+ P22 


Pl3 


+ P24 


P31 


+ P42 


P33 


+ P44 


Pll ~ 


f P33 


Pl2 


+ P34 


P21 - 


f P43 


P22 


■f P44 



(42) 



(43) 



(44) 



12 



The tensor product of two matrices pA and pB has the form quadratic in matrix elements of the matrix p 

PA ® PB 

The purity parameter of the two-spin state ((421) equals 



(Pll + P22)PB (/013 
(P31+P42)PS (P33 



P2i)PB 
Pii)PB 



E 



\pik\ 



i,k=l 

The purity parameter of the states ( ^3|) and (^ ) read 

PA = \pil + + |P13 + P2i\^ + |P31 + Pil? + 1/333 + Pii? 

and 



Mb = l/Oii + P33I 



|Pl2 + P34p + |P21 + P43p 



IP22 + P44| 



(45) 

(46) 
(47) 



One can calculate the trace defining the angle between two vectors corresponding to the density operators 
in the form 

Tr [ppA ® pb^ ^ 

[Pll + P22) Pll(pil + P33) + P12{P21 + P43) + P2l{pi2 + P34) + P22{P22 + Pn) 

+ (P31 + /O42) Pl3(Pll + P33) + PIa{P21 + P43) + P2z{Pl2 + P34) + P2a{P22 + P44) 

+ (pi3 + P24) P3l(Pll + P33) + P32(P21 + P43) + P4l(Pl2 + P34) + PA2{P22 + PAa) 

+ (P33 + P44) P33(Pll + P33) + P34(P21 + P43) + P43(Pl2 + P34) + P44(P22 + P44) 

(48) 



Having expressions (HQ), (45)^(|4q) one can calculate measure of entanglement for arbitrary density matrix 





/ 










1 





1 


1 





2 





1 


1 







I 








0; 



of two spins. For example, in the case of pure state 



P = 



one has 

, 1 1 

M = 1 ^-4=2 -"^ = 2 

This provides the maximum entanglement of the state (^), i.e. 

3 

For more general matrix of pure state of the form 



(49) 



(50) 



(51) 



/ 








M 







sc 








sc 







I 








/ 



(52) 



where c = cos ip and s = sin(p, one gets 



sin 2(^ 



The angle in ( |40| ) is determined by the angle Lp of (|53|). For Lp = 45°, one has 9 = 60' 



Cip = - sin 2lp 



(53) 
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12 Distance and entanglement of Gaussian squeezed states 



One can use the developed approach to study entanglement of two-mode squeezed Gaussian states. Wigner 
function of generic squeezed and correlated state p in n dimensions has the form {h = 1) [p8| 



W{Q)=Ne^p 



1 



o(Q)s~MQ) 



N = (detE) 



-1/2 



where Q = (pi,p2, ■ ■ • ,Pn, gi, 92, • • ■ , and 



Q = Q (Q) 



(54) 



(55) 



with (Q) being the parameters describing means of the quadature components. The 2nx2n matrix S 
describes variances and covariances of the quadrature components. We present the dispersion matrix in the 
form 



El Ei2 



^12 



A B 
B* C 



(56) 



We consider two subsystems with dimensions n\ and ^2, with n\ + n2 — n. Let us suppose that the system 
state has the parameters (Q) — (squeezed vacuum in the case of pure states) . The normalization constant 
N in (H) is determined by the matrix E due to the condition 



The purity parameter of the Gaussian state of the system equals |2S 



M = Trp2= J Vt^2(Q)iQ"^2-"(dctE)-i/2. 



Integrals (p7|) and (|5l) can be calculated using the formula for n-dimensional Gaussian integral 



-xax+bx 



dx 



•i/2 



: exp ( -ha 



vdet a 

The Wigner function of the subsystem state 1 pi, which is denoted as Wi{Qi), is given by the relation 

dQ2 



VFi(Qi 



(27r)"2 

and the Wigner function of the subsystem state 2 p2 is given by analogous relation 

M^2(Q2) = / W{Q)^. 
Both integrals are Gaussian ones. Due to this, one has 

Wi{Qi) = TViexp 

T^2(Q2) = A^2exp 

where 



-^(Qiar'Qi 



Q20'2 ^Q2 



Ni = (detcri)-^/^ 
N2 = (detCTa)"'/' 



(57) 
(58) 
(59) 
(60) 
(61) 

(62) 
(63) 



crfi = ^- SC^^B* a:^^ = C - B^A-^B. 
The purity parameters of the states of subsystem read 

Ail = 2-"i(detai)-i/2 p2 = 2""= (dot ^2)"^/^. 

In the case detE = (1/4)", fii = ^2- 



(64) 
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The normalization constants A'^i and N2 are functions of initial dispersion matrix S. The Wigner function 
of the state pi (E) P2 has the product form 



W'l2(Q) =M^l(Ql)W^2(Q2)- 



This form is also Gaussian one 

W^i2(Q) = iVi2exp 



-i(Q.-Q) 



N12 = N1N2 



0-1 

CT2 



(65) 



(66) 



To calculate introduced measure of entanglement, one has to calculate the fidelity t = Tr( ppi (81 P2 ) , which 
is expressed in terms of Wigner functions by the integral 



t 



dQ 

(2^ 



W{Q)Wi2iQ). 



(67) 



t= det (S + a] 

This trace determines the term with cos 9 in the expression for entanglement of the squeezed Gaussian state. 
Thus measure of entanglement of the squeezed Gaussian state reads 



The integral is Gaussian again with the dispersion parameters (cr + E) . So one has 

-1/2 



ec = 2-"(detS)-i/2 _^2-"(detcr)"i/2 - 2(det(i; + CT) 



-1/2 



It is determined by the quadrature dispersion matrix of the composite system, which is characteristic of 
Gaussians. If the pure state is squeezed but not correlated |2£, 27|, entanglement is absent. For entanglement, 
one needs the correlation of quadratures in the initial pure states.]] 



13 Purification of separable density matrix 

In this section, we consider the procedure of purification of mixed density matrix. Density matrix of a 
composite system is said to be simply separable if it has the form 

PAB ^ PA® PB (68) 

where pA = ^t:bPab and pB — Tr^p^^. Such matrix can be pure if and only if pA and ps are pure and 
hence paBi Pa and pB are projectors in the appropriate spaces. Such a pure state of the composite system 
is not entangled. More generally, a density matrix pab is said to be separable if 

PAB = ^Pn/0„A ® PnB ^Pn = ^ Pn > 0- (69) 

n n 

In this case, 



PA = ^^PnPnA and PB ^ 



PnPnB- 



Clearly if n > 2, the density matrix pab is not pure. For n — \ the matrix pab is not pure unless pA and 
Pb are one-dimensional projectors. Since pnA and pnB are density matrices, they are convex linear sums of 
projectors 

PnA^^Pnj^i-j and PuB^^quk^^k (70) 



^It is worthy noting tliat anothermeasure of entanglement based on the cross covariances of quadrature components of 
entangled modes was introduced in |35[. 
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where Pnj and Qnk are nonnegative numbers Pnj — X)fe Ink — 1- Thus pab is the convex sum of projectors 

PAB = PnjqnkUnjVnkuljvl^. (71) 

Note, these projectors are not all mutually orthogonal and pab is a mixture with weight PnjQnk for all n,j, k. 
The density matrices 

PAB {n, j, k) = UnjVnkuljvl^ (72) 

(there is no sum) are not mutually orthogonal pure state projectors, except for n ~ I. So, once the 
eigenvectors (eigenrays) are given, we could construct the density matrix of AB as a convex sum of projectors. 
For n = 1, what does pab contain that is not contained in p^gjps? In this case, it is a set of phase diferences 
between various eigenvectors that go to make up pab- There are rank(p/i) x rank(/9B) phases and hence one 
less phase differences. All this is very similar to the purification of mixed states. We write 

PA = Y1 Po^t PB = Y. 9fenf (73) 

j k 

and hence by hypothesis 

PAB= PA® PB=Ypjqk^f ®l^k- (74) 

jk 

So, for purification, we adopt the ansatz given earlier (we omit the Kronecker product symbol) 

1/2 nfnfn^^n^inf, 



= 2^ [PjikVi'qk') I =■ (75) 

jkj'k' /xr (n/nf n^sn^nf,n^sj 



Pab 



This pab is a pure matrix with probability weights pjQk for the projectors n^II^ and the suitable phase 
differences which number rank(p^) x rank(pB) — 1. It is essential to choose Hab such that Tr (n^^^n^IIj^) 
for all j, k. But the density matrix so constructed will not lead to pA and pb as partial traces. 
For a more general case of n > 2, we have the problem of rediagonalizing 

PA = Y^ PnjU^, PB^J2 l^k^nk ■ (76) 

n,j n,k 

Once this is done, we proceed as in the simply separate case (n=l) discussed above. 



14 Entanglement and straddling of fiducial projectors 

The fiducial projector in the purification protocol generates the relative phases between the (two or more) 

density matrices for pure states which have been 'superposed'. Such an operator is a 'phase correlator'. The 
Hermitian fiducial projector and the projectors which are being superposed are all Hermitian, yet relative 
phases are introduced. Given two or more Hermitian matrices one can generate the Bargman phase from their 
product. For this to occure, the operators cannot all commute (phase (ABC) 7^ 0). Simple examples may 
be provided by the Pauli matrices (or Dirac matrices), in our case, the phase difference between amplitudes 
is generated by the overlap of the fiducial matrix with the respective density matrices. Since we do a 
normalization in (12) only the phase of the overlap survives. There is a source of the phase interference 
introduced in our composition law (12). The question naturally arises — can we choose the fiducial projector 
Po so as to produce any set of phase differences? The answer is affirmative but not unique. 

The fiducial projector must straddle the pure states which are added, that is, it must have nonzero 
overlap with each of them. (It may or may not have overlap with other states.) Following up on this notion 
we find that the fiducial projector which restores a fully entangled pure state of a composite system straddles 
the eigenprojectors of the individual rays which are direct product of pure density matrices. (This is for 
separable systems, otherwise we get some direct product pure states and some fully entangled pure states.) 
This straddling implies that automatically the fiducial projector is a 'nonlocal' operator acting coherently 
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on the subsystems. Since it is also a projector, it follows that this projector is a fully entangled pure state. 
Only such an entangled projector can restore full entanglement. Fully entangled operators can be multiplied 
by each others or added together to obtain fully entangled operators, but they will not be projectors of rank 
one. This entanglement (and phase coherence) can be inherent in operators as well as in states. 



15 Conclusions 

We presented an intrinsic approach to different quantum phenomena which are entanglement and inter- 
ference. The approach is intrinsic because it points out the unique basis for both phenomena which is 
superposition principle of quantum states. But to use this superposition principle in generic case of mixed 
states, one needs the addition formula for density operators. The discussed measure of entanglement is 
intrinsically connected with the given state of a composite system because it is determined completely by 
partial traces of the state density operator and by the deviation of the density operator from the tensor 
product of the partial traces. Thus, because the entanglement is the property related to the state super- 
positions (expressed in terms of new addition rule of density operators with using a fiducial projector) the 
fiducial projector becomes a useful tool for treating both phenomena — interference and entanglement. 
To conclude, we point out new results of the paper. 

The nonlinear addition rule for impure density matrices, which results in pure density matrix given in ( p^ ) 
and (^4|), is a new purification procedure. The addition rule formula (|2^) for symbols of density operators 
of any kind (including Wigner distribution, tomograms, etc.) is another new result of our consideration. 

The notion of measure of entanglement of arbitrary order for bipartite and multipartite systems for an 
arbitrary observable given by (|4l| ) is a new aspect of entanglement suggested in our study. As a partial 
case, the measure introduced contains the description of measure of entanglement of density operator for 
bipartite system given in (^o|). The measure of entanglement is related directly to intrinsic properties of 
density operator of a composite system. 
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Appendix 1. Partial density matrices for pure state of composite 
system 

Let us consider the pure state of a composite system which has two subsystems A and B. The pure state is 
described by a vector of the form 

N M 

tp =^^^Cia(piXa (77) 
i—1 a— 1 

where is dimension of the subsystem A, M is dimension of the subsystem B, and the orthogonal vectors 
(fi {i — 1 . . . N) and Xa {a — 1 . . . M) form basis in Hilbert spaces of the subsystem states. 

The density operator of the pure state which corresponds to the decomposition ( |77| ) of the state vector 
has the form 

N M 

PAB^lpi^^^ ^ ^ CiaC*i3ipiip]Xc,xl- (78) 
ij = l Q,/3=l 

The density matrix of the A-subsystem state in the chosen basis has the matrix elements expressed in terms 
of decomposition coefficients 

M 

(p^)y =5]C,„q„. (79) 
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The density matrix of the iJ-subsystem state in the chosen basis has the matrix elements 

N 

{PB)a(} = ^ CiaC*fi. (80) 
i=l 

Both density matrices pA and pg are nonnegative Hermitian matrices and Tr p^ =Tr pg = 1. Let us calculate 
parameters 

/i(-^)=Tr(p^r A^i^)=Tr(pBr (81) 



for arbitrary integer n. 
One can easily see that 

In fact, 



(4f^=t4f^. (82) 



and 



JV M 

il,'i2,...*A/ — 1 CKl ,C^2,■■■tt7^ — 1 

N M 

f^ri = ^ ] y ^ ^iiai ^iiQ2 ^«2a2 '-^igas ' ' ■ C'i„_ia„„iC'j^_j^jj^C'i^Q,„ C'i^Qi • (84) 

■il ,22 ,. . . in — 1 Ctl ,Ct2 , ■ ■ -CKti — 1 

The terms without star are the same in both expressions ( p3| ) and (p^). These terms are invariant with 
respect to arbitrary permutations 

1,2, . . , ,n ^ Sl,S2, ■ • • ,Sn- 



The terms with star look differently in ( |83| ) and (84), but since both sums ( |83D and (84) are invariant with 
respect to arbitrary permutations, let us make the particular permutation 

1, 2, . , . , n — 1, 71 — > 2, 3, , . , n, 1 



in sum (g4|) . The terms without star are invariant and the terms with star in ( p4D after permutation coincide 
with the terms with star in (|8^) . This proves equality (|8^) which means that eigenvalues and rank of the 
matrices pA and ps are the same. It is clear that the proof can be extended to multipartite composite system 
AB . . .C. Thus we get the following result. Given a pure state of a multipartite quantum system pab...c- 
Then the eigenvalues and ranks of the density matrices pA, Pb, ■ ■ ■ Pc are equal. 
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